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ABSTRACT 

We derive the fully relativistic Ohm's law for an electron-positron plasma. The absence 
of non-resistive terms in Ohm's law and the natural substitution of the 4-velocity for the 
velocity flux in the relativistic bulk plasma equations do not require the field gradient 
length scale to be much larger than the lepton inertial lengths, or the existence of a frame 
in which the distribution functions are isotropic. 
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For a plasma, Ohm's law describes the relation between the induced current and the 
plasma electric field. For an ion-electron plasma, the field depends on resistive, inertial, 
and Hall effect contributions. The result is usually derived for the non-relativistic limit in 
the Boltzmann picture 1 . Ohm's law plays a direct role in the magnetic induction equation 
used in the description of bulk plasma dynamics. 

Relativistic plasma models have been effective in explaining the observations of rela- 
tivistic astrophysical jets and winds 2 . Such models have generally employed the relativistic 
continuity equation as a vanishing 4-divergence of the bulk 4-velocity, and Ohm's law as a 
simple bulk 4- vector generalization of the non-relativistic equation 2 ' 3 . 

We shall see that for a two-component relativistic plasma composed of different mass 
particles, the natural use of these magnetohydrodynamic (MHD) forms for the continuity 
equation and Ohm's law requires the existence of a reference frame in which both distri- 
bution functions are isotropic in momentum. The constraint results from the non-linear 
relation between momenta and velocities in the relativistic regime. 

This requirement is non-trivial because distribution function isotropy also requires 
the plasma under study to be microinstability saturated; otherwise microinstabilities could 
grow because of distribution function anisotropy. Yet, evidence for the presence of anisotropies| 
and microinstabilities in relativistic winds has come from observations of the Crab Nebula 4 . 
In particular, the Weibel instability has been suggested to explain the presence of wisps 
downstream from relativistic shocks 5 . Anisotropies downstream from relativistic shocks 
are likely present in jets as well. 6 
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Depending on density and temperature conditions, anisotropies in relativistic plasmas 
may also arise from anisotropic radiation onto a plasma, such as in coronae of stars or in 
models active galactic nuclei (AGN) for which a pair atmosphere forms. 7 The impinging 
of winds and jets onto ambient media also produces anisotropies 8 . 

In this paper we show that in contrast to an ion-electron plasma: 1) only a resistive 
contribution to Ohm's law for a relativistic e + — e~ plasma is relevant under quite general 
conditions, and 2) that anisotropy in the distribution functions need not affect the form 
of the unperturbed relativistic bulk equations for the pair plasma. Thus, relativistic bulk 
dynamics for pair plasmas which exhibit evidence for microinstabilities are appropriately 
described by the relativistic MHD formalism, whereas ion-electron plasmas which exhibit 
such instabilities are not. 

The Boltzmann equation is given by 



df/dt + v*dj + F*d p J = [A t f] colL + [A t f] 

radiat. creation annihilation (1) 



where v l is the particle velocity, p l is the particle momentum, t is the time, F l is the 
electromagnetic force (since we ignore gravity), / = /(x, p,t) is the scalar distribution 
function, and the terms on the right are schematic. We shall assume that collisional losses 
dominate synchrotron radiation losses, which is acceptable for 10 



where 7 is the particle Lorentz factor, is the pitch angle, e is the positron charge, and 
B is the magnetic field measured in Gauss. 




(2) 
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Define the plasma quantities: number density 

n± = J /±(x,p,t)d 3 p, (3) 

and velocity flux density 

A = J v* ± f ± (x,p,t)d 3 p = n ± (v* ± ). (4) 

Adopting the signature (+, — , — , — ), we then have the flux density 4-vector 0± = (cn±, 0±),| 
and current density 4-vector j M = e((p fJ ' + — </> M _). The energy and momentum densities are 
the (00) and (iO) components of the symmetric kinetic tensor given by 



K ik 



J v* ± p k ± f ± d s p = n ± (v* ± p k ± ) (5) 



and 

K^ = (6 ± ,cir ± ), (6) 

where U± is the momentum density of either positrons or electrons. 

Quantities without the ± shall refer to the plasma as a whole- the sum of the compo- 
nent contributions. We require the existence of a proper frame moving with 4-velocity 
in which the charge density, the velocity flux density, and the momentum density vanish. 
We denote quantities in this frame with a superscript *. For an e + — e~ plasma this means 

j* = e(n* + -n*_) = 0, (7a) 

r = 0, (76) 

and 

IT* = m e (n* /2)« 7+ <>* + <7-<0*) = 0. (7c) 
4 



Since K^ v has 10 independent components, we can write 



X* ■— P* 4- A* 



(8) 



where the symmetric tensors P^ v and satisfy 



P* 1 =P5 tJ , P * = 0, P o = e*, 



(9) 



and 




0. 



(10) 



In (8) — (10), P is the scalar pressure and A*- has 5 independent components that measure 
anisotropy. 

Finally, define the 4-vector by 



Thus from (76), H l * pair vanishes. But H° pair = by definition, so H pair * = 0. Since H pair 
is a 4-vector, the vanishing of H^* pair implies that H^ air = in all frames. Note that H^* 
measures the heat flux density per unit mass in the proper frame, so we shall call the 
heat flux density 4-vector. 

The procedure used to derive Ohm's law for a relativistic pair plasma is as follows: 
(i) First we obtain a "resistive" type collision term appropriate for a nearly collisionless 
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H" = [2/(m+ + m_)] - (m+(f>l + m_<^)]. 



(11) 



where p* is the proper frame rest mass density. For a pair plasma this becomes 



C = (lK)]^-me(f)]. 



(12) 



relativistic pair plasma, (ii) Second, we relate this to the current density, (iii) Finally, a 
subtraction of the momentum density equations for the positrons and electrons yields the 
desired result. 

We can find the value of the momentum density in any frame by Lorentz transforming 
the stress-energy tensor. The result is 

clT = c 2 m e (n+(u\) +n_(w l _)) = 7 j>(P + e^V'c' 1 + A°\ (13) 

where u l is a spatial component of the particle 4-velocity, V 1 is a component of the bulk 
3-velocity, 7y is the bulk Lorentz factor, and the anisotropy term is given by 

A 0l = A* j U j + U k U l A* kl uy( 7v + 1). (14) 

Note that since the proper frame 4- vector A^ V *U* = 0, we know that it is zero in all 
frames. Thus A i0 = A ik U k and (13) can be written 

clT = c 2 m e (n + (u\} + n_(u 1 _}) = <y v (P + e*)U l + A lj U 3 . (15) 

Inverting (15) we obtain 

lv V'c- 1 =U'= [c 2 m e (n+(u\) +n_ <«!_)) - ^-]/[( 7 y)(P + e*)]. (16) 

The average 4-momentum gains from collisions in the proper frame are given by 

Af/t* = -A P r = (m e /2)(«>* - (</>*). (17) 

The approximate proper frame pair plasma collision term is then 

(nV c */4)(Af/n = -(nV c 74)(A/_*) = = -P» + \ (18) 
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where v* is the proper frame collision frequency. 

Equations (12) and (16) give 

m e c 3 (n + (u\_) + n_(w l _)) - cA lJ U 3 = 7v (P + e*)n* _1 (0!_ + </> l + ). (19) 
Now in the electron frame, we have 

c 2 n; (_) - [cA^Uj]^ = 1 ^\P + e*)n*- 1 <f)\_ { ~\ ( 20 ) 
and in the positron frame 

c 2 ir_ (+) - [cA^Uj]^ =7{, +) (P + e*)n*"V- (+) - (21) 
Transforming (20) and (21) to the proper frame, subtracting, and using (7) and (18) gives 

z/ c *(IT; - ITI) = 2P*1 = n W, (22) 
where the effective resistivity r\ r is given by 

Vr = K^er'c-'up+en^er^fh: + i)-y*.v(en* C )- 2 + 7c * + 1] 

^VCen*)" 1 - 2i/*- 1 7 c * 2 (7c* + 1) _ V«0 _ V^*}» (23) 
and 7 * is the Lorentz factor corresponding to the velocity 

(n*/2)-y ; = -(n*/2)-y *_ = /'/(en*). (24) 

We have used the 4- vector indices in r\ r since jo* =0. 

Equation (22) suggests that we subtract the momentum equations for the electrons 
and the positrons to obtain Ohm's law. This is standard in the non-relativistic case, but 
is only fruitful in the relativistic case because H^ air = 0. 

7 



Under the conditions of (7), conservation of energy and momentum momentum implies 
that the contributions of the annilihation and creation terms of (1) to the the energy and 
momentum equation cancel for each plasma component in the proper frame. Thus, by 
covariance, the contributions to the 1st moments in all frames vanish, and we do not 
consider them further. 

The zth component of the proper frame relativistic energy mometum tensor for positrons,! 
as obtained from the 1st moment of the Boltzmann equation is given by 

<9 IT + * = -d k Kf - e(n*/2)E** - [e(v+/c> x B]** + (25) 

Subtracting the analagous equation for electrons, and using (8), (18), and (22) we get 

E » = Vr r - m e c(2ep*)- 1 a fc (Aj* - A k f) - <9 [(W^).f *]• (26) 

The requirements (7a) and (76) which led to (26), are assumptions about the Oth and 
1st moments of the Boltzmann equation. If we further assume that 

Af = A l l\ (27) 

the second term on the right hand side of (26) would vanish. Then, in the steady state, 
we would be left with 

E H = rjrf*. (28) 



Note that the dependence on the current in (28), results from (22). The latter follows 
for example, even in the Fokker-Planck approximation if, as in our case, the charge density 
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vanishes in the proper frame. In addition, note that none of the assumptions that led 
to (28) require the existence of a reference frame in which the distribution functions of 
positrons and electrons are isotropic. 

In a general frame, (28) becomes 

F^U„ = Vr (j»+fU T Un, (29) 

where we have added a convection term. Equation (29) resembles the magnetofluid relation 
11 . Here however, the effective resistivity r\ r depends on a scalar function of the current 
and momentum density of the plasma components. This dependence is eliminated when 
we choose the coordinate axes such that the proper frame velocity lies on a principal axis. 
In addition, 7* factors out of (23) and we have 

Vr -> r,^) = v* c c- 2 (n*e)- 2 {2(P + e*) - e*} = z/ c *c" 2 \n* e)" 2 ! (2P + e*), (30) 

where the superscript (pr) indicates that the current flows along a principal axis. 

Note that unlike the Ohm's law for an ion-electron plasma, there are no Hall effect 
or pressure contributions to (29). The vanishing of the bulk Hall effect term is the result 
of an effectively zero net particle gyration frequency; the sum of the electron and positron 
contributions vanish due to opposite streaming of positrons and electrons around the field 
lines. The vanishing of the pressure term results because our assumptions have appealed 
to the mass symmetry of the problem, eliminating diffusion. The Hall effect and pressure 
terms are negligible for an ion-electron plasma only when Xf g » Aj, and Xf g » \j3, 
respectively, where A f g is the length scale of the field gradients, A; is the ion-inertial length, 
and P is the ratio of thermal to magnetic pressure. 
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The vanishing of the heat flux, which led to (22), is also the condition which allows 
substitution of the bulk 4- velocity components for the bulk 4-flux components in the plasma 
continuity equation. The point is that the usual integration of (1) over momentum gives 
dfj,^ = 0, so substitution of the 4-velocity for the 4-flux 4> requires = 0. If ^ 0, then 
this substitution would produce an inhomogeneous equation. Because m+ >> m-, for a 
relativistic ion-electron plasma, H^* only naturally vanishes if both distribution functions 
are isotropic in the proper frame. This is true whether or not the proton component has 
a relativistic temperature. 

Note that vanishes for any two component plasma in the non-relativistic limit. 
To see this note that 

H** oc -m+n+(v i + )* + -m_ra* = m+n+d+v^. - v\)* + m_ra* (7_vl (31) 

where the second equality follows since the proper frame momentum density vanishes by 
definition. Thus when n* + = n*_, (31) vanishes when 7 + = 7_ = 1. Therefore vanishes 
by the argument that follows equation (12). 

Finally, note that our proper frame result is not the full nonrelativistic limit, because 
of the relativistic temperature. In the non-relativistic limit, we have the additional result 
that the pressure drops out of (23) and that e* = n*m e c 2 so 

rj(P r ) ^ v = (v c n*m e c 2 )/(n*ec) 2 (32) 

= m e u c / (n*e 2 ) = 2fi pair u c / (n*e 2 ), (33) 

where r\ is the non-relativistic resistivity, and fx pa i r = m e /2 is the reduced mass for the pair 
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plasma. Equation (33) is the usual non-relativistic result. In the case of a non-relativistic 
ion-electron plasma, [i va i r is replaced by Hi e = mim e / (rrii + m e ) ~ m e . 

We would like to thank A. Loeb for discussion. 
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